Emission characteristics of quantum dots in planar microcavities 
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The emission properties of single quantum dots in planar microcavities are studied experimentally 
and theoretically. Fivefold Enhanced spontaneous emission outside the microcavity is found for dots 
in resonance with the cavity mode, relative to detuned dots, while their radiative lifetime is only 
marginally decreased. Using high power excitation we obtain the in-plane cavity dispersion. Near 
field images of the emission show spatial distributions of several microns for resonant dots, which 
decrease in size with the detuning from resonance. These experimental findings are explained using 
a simple and intuitive model. 
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I. INTRODUCTION 

The potential applications of semiconductor quantum 
dots (QDs) as quantum light emitters have generated 
considerable research efforts in recent years i Single pho- 
ton sources using quantum dots are an important ingre- 
dient in quantum information applications such as quan- 
tum cryptography and teleportation,^ and may be fur- 
ther employed to implement efScient linear optics quan- 
tum computation"^ with the attainment of on-demand 
indistinguishable single-photon pulses."* Single QDs can 
be excited either with continuous wave (cw), where 
they emit antibunched photons obeying a sub-Poissonian 
statistics or with pulsed light, where a single photon is 
emitted per pulsei It is the radiative recombination of 
electron-hole (e-h) pairs (or excitons) confined in single 
QDs that gives rise to the emission of single photons. 

QD based devices have several advantages as single 
photon emitters. These include relatively large oscillator 
strengths, narrow spectral linewidths, and compatibility 
with mature semiconductor technologies. A major hin- 
drance in their usefulness, though, is the low extraction 
efficiency of the emitted photons due to the high refrac- 
tion index of the host semiconductor. Typically, only 
about two percent of the photons will be emitted from a 
GaAs based device, while the rest will be lost due to to- 
tal internal reflection. The general approach to overcome 
this problem has been to place the QDs at the antin- 
ode of a microcavity whose dimensions are comparable 
to the wavelength of the emitted photons.'^ In these cav- 
ities the number of allowed optical modes is reduced and 
at the same time, the in-cavity intensity of the allowed 
modes increases. As a result, the spontaneous recombi- 
nation rate of excitons within these resonating modes is 
increased, a phenomenon known as the Purcell effect &i 

Several approaches to realize 3D photon confinement 
have been studied, including whispering-gallery modes 
in microdiskSfS defect modes in 2D photonic crystals)^ 
and lateral patterning of planar dielectric cavities using 



electron-beam lithography and etching J^IS The latter has 
proved particularly useful, demonstrating efficient pho- 
ton emission in the weak coupling regime^ and more re- 
cently strong coupling, paving the pathway for solid state 
realization of coherent control schemes. In the strong 
coupling regime, the interaction between the cavity mode 
and the emitter is larger than their combined decay rates, 
and the irreversible spontaneous emission from the QD 
exciton is replaced by a coherent exchange of energy be- 
tween the exciton and the cavity mode. Among the chal- 
lenges in the implementation of strong coupling in struc- 
tured cavities, one notes the reduction in the cavity's 
Q-factor due to the lateral patterning, and the lack of 
control in the in-plane positioning of the QD with respect 
to the lateral microstructure, resulting in a reduction in 
the oscillator strength of the confined excitoniiStii 

In view of the technological difficulties associated with 
systems of QDs in 3D cavities, it seems important to 
study a system of QDs embedded in a planar microcav- 
ity, which is easier to fabricate. These cavities consist of 
distributed Bragg reflectors (DBRs), which are typically 
stacks of alternating quarter-wavelength-thick layers of 
GaAs and AlAs, separated by a spacer layer of GaAs. 
Planar cavities can attain high Q factors due to the excel- 
lent control in their fabrication and embedding the QDs 
at the field's antinode is straightforward. These struc- 
tures support an in-plane continuum of modes, and the 
interaction with QD excitons is therefore not expected 
to be in the strong coupling regime Nevertheless, the 
spontaneous emission and the light extraction efficiency 
can still be enhanced considerably.^'^ For example, a rel- 
atively weak planar cavity has been utilized recently to 
produce a ten-fold increase in the efficiency of light col- 
lection of GaAs based LED structuresii^ 

Unlike systems of quantum wells embedded in planar 
cavities, where polariton effects prevail and have been 
extensively studied (see e.g.fi^*iSiiLi&), QDs in planar 
cavities have received little attentioniiSi^" Among the 
few experimental studies of this system, angle-resolved 
photoluminescence spectroscopy was used to detect TE- 
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TM mode splitting (not to be confused with Rabi split- 
ting, which is discernable only in the strong coupling 
regime) ^ 

In the present work we study experimentally and theo- 
retically single QD emitters in a planar microcavity. Us- 
ing a partially covered layer of InGaAs self assembled 
QDs placed at the center of a A microcavity, we have mea- 
sured various characteristics of the spontaneous emission 
dynamics of the system such as spatial and angular distri- 
butions of the radiation field and its temporal decay. Our 
results are compared with a relatively simple theoretical 
analysis, which follows the treatment of Sugawara^^ and 
gives a simple and intuitive picture of the spontaneous 
emission properties of recombining excitons confined in 
QDs embedded in a planar microcavity. This approach 
accounts reasonably well for our experimental findings, 
including emission rates, radiative lifetimes, and angu- 
lar and spatial distributions of the emission. In all of 
these measurements, a crucial parameter is found to be 
the detuning of the emitter energy from that of the cavity 
resonant mode. Since we do not have an exact knowledge 
of the QDs size, shape and composition, detuning from 
resonance is modeled in our approach by either varying 
the QD size or by varying its material bandgap, while 
keeping the optical properties of the microcavity fixed. 

The paper is organized as follows. In section^]we pro- 
vide the details of our experimental system. Section lllll 
gives the theoretical model including the exciton wave- 
functions and their interaction with the cavity field. In 
section IIVI we give the experimental measurements and 
compare them with model calculations. A summary is 
given in section and some details pertaining to the 
calculation of the excitonic wavefunction are given in the 
appendix. 



II. EXPERIMENTAL SETUP 

The samples were grown by molecular beam epitaxy 
on a (100) oriented GaAs substrate. One layer of strain- 
induced InGaAs QDs was deposited in the center of a 
one wavelength GaAs spacer layer. The height and com- 
position of the QDs were controlled by partially cover- 
ing the InAs QDs by 30A thick layer of GaAs and by 
subsequent two minutes growth interruptions^ to allow 
diffusion of In (Ga) atoms from (into) the strained is- 
lands. The growth resulted in InxGai_xAs QDs whose es- 
timated radius and composition are 150 A < Rqd < 250 A 
and 0.47< x < 0.53, respectively, and height of 30A. The 
GaAs spacer layer was grown to a width whose resonance 
mode matches the QDs emission due to ground state e-h 
pair recombinations (1 A cavity). The microcavity was 
formed by 25 and 11 period stacks of alternating quarter 
wavelength layers of AlAs and GaAs below and above 
the spacer layer, respectively. The samples were not pat- 
terned or processed laterally to prevent obscuring the 
emission and its polarization. 

For the optical measurements we used a diffraction lim- 



ited low temperature confocal optical microscope^i*^. 
The sample was cooled by a copper braid attached to 
the cold finger of a He- flow cryostat. The sample mount 
was accurately manipulated in three directions using 
computer-controlled motors. A X60 in-situ microscope 
objective was used in order to focus cw or pulsed laser 
light at normal incidence on the sample. The emitted 
light was collected by the same microscope objective. 
The collected light was spatially filtered, dispersed by a 
0.22 m monochromator and detected by a nitrogen-cooled 
CCD array detector. The system provides diffraction- 
limited spatial resolution, both in the excitation and the 
detection channels. For the time-resolved spectroscopy, 
the dispersed light from the monochromator was focused 
onto a small, thermoelectrically cooled, single channel 
avalanche silicon photodiode. The signal from the photo- 
diode was analyzed using conventional photon counting 
electronics. The photodiode dark count rate was 100 
sec""'^, and the system temporal resolution was approxi- 
mately 250 psec. 

III. THEORETICAL FRAMEWORK 

The following approach assumes a single Is exciton 
confined in a QD interacting with a discrete cavity mode. 
No charged exciton or biexciton states were considered, 
although their presence was confirmed in the measure- 
ments (see section FlV A|) . Our main interest in the cur- 
rent study is in the cavity effects on the excitons emission 
dynamics. The assumption of discrete cavity modes is 
justified for high finesse cavities, (we have verified that 
only the basic A mode contributes appreciably to the in- 
teraction with the exciton) . Although more complete de- 
scription of the electromagnetic density of states in the 
cavity is availableii, the main features of the exciton- 
photon coupling are well demonstrated within the frame- 
work of our simplified model. 

A. Model wavefunction for quantum dot excitons 

We start by solving for the excitonic wave function in 
a quantum dot, considering finite potential barrier in the 
growth (z) direction and parabolic potential in the lateral 
direction. 

The effective mass Hamiltonian is given by 



<^\/{Pe - PhY + {Ze - ZhY 

Ve{pe,Ze)+Vh{ph.Zh), (1) 

where = (p^, Zi) are the electron or hole in-plane and 
z coordinates, and a = ml/m)^, a^z = mi/ml^, auz = 
mi/m'l^ are the appropriate effective mass ratios in the 
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plane and z directions, and e is the background dielectric 
constant screening the Coulomb interaction. In Eq. 
we have neglected the differences in the conduction- and 
valence-band masses and in the dielectric constant value 
of the two semiconductors that comprise the QD and its 
host. In the case of strong z direction confinement, the 
electron and hole potentials can be approximated as 



(2) 



thus decoupling the problem into lateral and z parts, and 
enabling us to write the envelope wave function for the 
exciton in the form 

(for infinite z direction potential, this approximation is 
exact). Transforming to center of mass (CM) and relative 
in-plane coordinates: 



R = — j— ; Pe- Ph, 

1 + a 



(4) 



we separate the CM-motion and relative-motion wave 
functions, \I'x(R, r) = ■;/'(R-)0(r), by Choosing parabohc 
potentials for the lateral confinement of both electron 
and holeS^ 



1 



ViiPi) = -jjn'i^ Pzi i = e,h. 



(5) 



Employing natural units of length and energy, namely, 
bulk effective Bohr radius (a^ = h'^e/mle^) and Rydberg 
(Ry — fi? / 2m\a?g) ^ the Hamiltonian in Eq. takes the 
form 



H — Hji + Hr, 



Ha = V 



1 + <7 ■ « ^ i\\^a) 

Hr = - [(1 + + (Jezdl + CJh^dl^ 



(6) 

-^V^{z^)^Vh{zh) 



where ^ is a dimensionless lateral localization parameter 
related to the confining potential by 



as 



(7) 



and Afx = rUe + m|[ is the in-plane CM exciton mass. 

We identify the in-plane CM motion part with a 2D 
harmonic oscillator Hamiltonian, which is solved analyt- 
ically resulting in wave functions '0(R) which are given 
by the associated Laguerre polynomials^'"^. In what fol- 
lows, we will be interested only in the CM-motion ground 
state which is a 2D Gaussian function. By equating the 
areal size of the CM wave functions with the physical 
QD radius: \ J ij{Ti)(fR\^ ^ 2n^^ 
potentials |SJ| with the QD size. 



7ri?|), we relate the 



For the in-plane relative coordinates we employ a trial 
wave functions^ 



(r) = AA,e-('-/'')° 



(8) 



where Afr is a normalization constant and 77, a are vari- 
ational parameters, a has values between 1 (no lateral 
confinement - very large QDs) and 2 (strong lateral con- 
finement - very small QDs). Considering finite height 
potentials in the z direction, Vi{zi) — V{'Q{\zi\ — L/2) 
where i = e,h and L is the QD height, we have for the z 
wave functions^SiSi 



Xi{z^) =7V;, X 



cos{kiZ,) 



ki 



-K.(lz.|-L/2) 



\z,\ < L/2 

\z,\ > L/2 (9) 



where i — e,h, Mz^ = (^/2+ l/^i) is a normaliza- 
tion constant. 



2m\ 



{Vf - E^,), (10) 




and E^Q is determined by the ground state quantization 
condition: 



(11) 



The variational parameters in Eq. (jSJ are calculated by 
maximizing the magnitude of the ground state exciton 
binding energy Eg = E'^„ + E^^ + E'^^ - E,, where 
= fiu! = ^2 (1+cr) ^^'^ ground state energy of the 
CM motion and E^ is the expectation value of the relative 
motion part of the Hamiltonian in © (see Appendix). 

The InGaAs self-assembled QDs that we are con- 
sidering have an estimated height of L = 30A. The 
conduction- and valence-band in-plane and normal to 
the plane masses were calculated for this quantum well 
width, assuming 50% InAs contents, using an eight 
band Kane model, and taking into account both non- 
parabolicity and lattice-mismatch strain effects.^'* We 
find for the electron: ml = 0.0630, ml = 0.0897m, and 
for the heavy hole: m| — 0.1573m, m^ = 0.3406m. 
For the potential barriers we find: V/ = 441.9meV 
and Vf^ = 110.6meV. Using these values we minimize 
Eq. (jA-ip with respect to the two variational parameters 
a and rj, and find the binding energy of the ground state 
heavy-hole exciton as a function of the QD radius shown 
in figure ^ 



B. Exciton-photon interaction in a planar 
microcavity 

In the following we assume weak coupling between 
the QD exciton and the cavity modes. This is justi- 
fied a-posteriori by the calculated spontaneous emission 
linewidths, which are always much smaller than that of 
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given (using the Coulomb gauge) by: 




Dot Radius [A] 

FIG. 1; (a) Exciton binding energy, (b) Variational parame- 
ter a, (c) Variational parameter 77, and (d) Detuning of the 
exciton line with respect to the cavity mode, as functions of 
the QD radius (resonance was assumed for Rd ~ 200A). 



2Trhv 



cos(/fc^z)e|| [e*ii-Pak + e"'"^" '^aHsa) 



2'Khv 



k " k 



(13b) 



where v — c/rics is the velocity of light in the cavity 
medium, Vc is the cavity quantization volume, and aj^ 
(dk) is the creation (annihilation) operator of the k field 
mode. In Eqs. H13|l kz satisfies the resonance condition 



kz = 



2rmT 



(14) 



where m takes any integer value and Lf, is the cavity 
width. Here we have considered only even modes which 
are the only ones which couple with the exciton ground 
stateji^ Assuming a system of a single exciton and a 
single photon, appropriate for low exciton densities, the 
time-dependent state of the system is given by 



IvE'(t)) =6(t)e-'"==*|$„0> + 



k,cr 



Ck, 



.(t)e-'"'=*|$g,W) 



(15) 

where to-x (tOk) is the exciton (electromagnetic field) res- 
onant frequency, and a runs over the two light polariza- 
tions. Taking the system to be initially with one exciton 
and the radiation field in the vacuum state [6(0) — 1, 
Ck,<T(0) = 0] we write the equations of motion for b{t) 
and Ck,CT(i), using the rotating wave approximation: 



the cavity. The cavity photon escape rate was deduced 
from the measured reflectivity of the cavity. We have 
put the cavity structure into a linear dispersion model 
and calculated it's reflectivity using transfer matrix for- 
malism. The resulting cavity mode linewidth was close 
to the measured value. 

The interaction between an exciton and the electro- 
magnetic field is given by 



int = -— Ep-A£ 

k,CT 



where AJ is the vector potential associated with the 
fcth mode with polarization a and we have neglected the 
quadratic term in A. In a planar cavity, with z in the 
growth direction, the Maxwell equations with the appro- 
priate boundary conditions yield two solutions for the 
electromagnetic field.-^ These correspond to the electric 
vector being normal to (TE) or in (TM) the plane de- 
fined by the wave vector k = (k||,A:z) and z, and are 



kt) - -i5].9k,.e'("^--^)*ck,.(i) - fb{t) (16a) 

k,(T 

cwAt) = -i.gk,,e-'("^-"-)*&(i) - YCk,.(0- (16b) 

In Eqs. Hl()|l 7x is the exciton broadening due to all pro- 
cesses other than spontaneous emission (phonon scatter- 
ing, nonradiative recombination) , and 7c is the decay rate 
of the cavity photon mode due to mirror losses, which is 
inversely proportional to the cavity quality factor. The 
exciton-photon coupling constant appearing in Eqs. H16|l 
is given by 

hgk.a = -($x, 0|— p • A£|$g, Ik,.). (17) 
mc 

The coupling constant is related to the total oscillator 
strength per unit area /k through 

2^ \gk.a\ = n-/k- (18) 

(T=TE,TM °" 

Using the dipole approximation (kg w — k;i), considering 
both polarizations and neglecting inter-subband mixing. 
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/k is calculated to be 



IV. RESULTS 



/k = lV'(fc|l)P 

where 

/ll = 

fi- = 



\mi 



(19) 



\Pc.F,,?\{Xe\Xh)? (20a) 

l^cvFej'KXelX/.)!'. (20b) 



In Eqs. (|19I20|I V'(fc[| ) is the in-plane Fourier transformed 
exciton CM wave function, 0(0) is the in-plane exciton 
relative motion wave function at zero, Pcv is the bulk 
transition matrix element, and i^gy (^e^) is the parallel 
(perpendicular) polarization factor in quantum wellsii^ 

Integrating Eq. H16b|) and substituting the resulting 
Ck,(T(i) in Eq- (|16a|l we have the integro-differential equa- 
tion 



h{t) 



1, ^ -^0 



-i(wk-Wx)-7o/2](t-r) 



■■-b{t). 



(21) 



In the weak coupling regime, which is clearly our case, 
7c ^ 5 and one can solve Eq. (|21|l to a good approxima- 
tion by taking h{t) out of the integral. For i 3> 7^7^ the 
result is 



where 



7SE = 



27r^'|gk,o-P>C(tJk -t^x, 7c) 

k,(T 



(22) 



V'— /:(wk - Wx,7c), 

^ Wk 



and 



£(wk 



,7c) 



1 



2^ (c^k - u^Y + (7e/2)2 



(23) 



(24) 



is the normalized Lorentzian cavity-mode broadening. 
The prime over the sum in Eq. (|23|l indicates taking the 
lowest cavity mode, — 2t: / L^, and summing over the 
in-plane photon wave vectors, k|| gx, where gx = ^x/v 
is the maximum exciton in-plane wave vector that can 
still couple to the cavity mode resulting in radiative re- 
combination of the exciton. The solution of Eq. (|22(l gives 
an exponential irreversible spontaneous emission of the 
exciton due to its coupling with the cavity field modes, 
representing Fermi's golden rule. As long as the calcu- 
lated 7sE is much smaller than the cavity broadening 7c, 
the weak coupling approximation holds; otherwise, one 
must solve the general equation H21|l . 



We are interested in particular in the influence of the 
cavity on the exciton's spontaneous emission, radiative 
lifetime, and emission distribution. 



A. Exciton Spontaneous Emission 

Figure [21 shows the measured reflectivity of the cav- 
ity and the integrated photoluniinescence intensity, using 
relatively high power, constant excitation. Each point in 
the flgure results from emission from a distinct QD. The 
numerical aperture of our microscope objective in the PL 
measurements was 0.8, resulting in a collection angle of 
~ 13.2° inside the sample. The flgure shows a sharp dip 
in the reflectivity corresponding to the resonant cavity 
mode. The FWHM is 1.94meV giving a Q-factor of 670. 
The emission intensity shows a pronounced enhancement 
for QDs in resonance with the cavity mode. A rapid de- 
crease in the emission intensity to a value roughly five 
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FIG. 2: (color online) Measured reflectivity (blue line) and 
integrated photoluminescence intensity (red points) of QDs 
at the antinode of a planar A cavity. 

In order to discuss these data, we relate the detuned 
exciton line to the A cavity width bjiiS 



A. 



-c 



277 

9x 



(25) 



where corresponds to the detuning of the exciton line 
from resonance. We assume the cavity resonant mode 
matches the emission from the heavy hole Is exciton 
of a QD with a radius Rd — 200A, for which ^ = 1. 
The corresponding resonance energy is calculated using 
Eq. (foil to be E'^'' = Ea = 1.2756eV. In order to model 
microscopically the emission from quantum dots that are 
detuned from this resonance, we represent the detun- 
ing by either changing the QD size or by changing the 
QD material bandgap (via the dot composition). The 
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dependence of the detuning on the QD radius is given 
in Fig. Considering only the basic cavity mode ex- 
pressed as kz — Qx/C, Eqs. (|^ and (fH^ are used to 
give 

C(v^q^ + ql/e-u:x,lc) ■ (26) 

In deriving Eq. we assumed that the QD resides in 
the center of the cavity, where the cavity mode takes its 
maximum value, and considered only the heavy-hole ex- 
citon whose perpendicular transition matrix element van- 
ishes, f± — 0. To account for the experimental collection 
angle, the upper limit of the integration in Eq. I|26(l is 
replaced with (/o = sin 0, where 9 is measured from the 
z-axis. 

In the calculation we used ricff = 3.5, e = 13.9 result- 
ing in QW exciton Bohr radius of = 116. TA. The free 
(QW) exciton oscillator strength per unit area was calcu- 
lated for a 30A Ino.53Gao.47As QW and was found to be 
/|| = 7.1 • lO-^A-2,^^ This value increases as the QD 
size reduces, due to the shrinkage of the relative motion 
wave function (see Eas. 1201) . We took an additional fac- 
tor of two for the oscillator strength to account for the 
two spin configurations. 

In Figure 13 the calculated exciton spontaneous emis- 
sion rates are shown together with the measured PL 
values as functions of exciton energy. The solid (blue) 
line corresponds to a calculation where detuning of the 
exciton energy was obtained by varying the QD radius 
whereas for the dashed (green) line the QD size was fixed 
to i?D = 2OOA and its bandgap was varied. In order to 
put the measured and calculated values at the same fig- 
ure, a scaling factor of 3000 was taken between the right 
(photon counts) and left (emission rate) axes. This fac- 
tor corresponds to the photon extraction efficiency in the 
experiment, which was estimated from pulsed excitation 
measurements!!^ 

The calculated results obtained by varying the dot size 
agree qualitatively with the experimental data, including 
the rapid fall off in the PL intensity away from resonance. 
The range of 90A < Rd < 300A that accounts for a de- 
tuning range of ~ 20meV is consistent with variations in 
typical self-assembled InGaAs QDs. The two model cal- 
culations in Fig.|21show marked differences. In particular, 
the emission enhancement shown for negative detuning 
in the case of variable QD size is missing in the fixed QD 
size calculation. The latter case shows a rapid inhibition 
of the emission as the QD becomes negatively detuned 
out of the cavity mode linewidth, which is clearly a cav- 
ity effect. In the case of variable i?D the cavity effect 
competes with the localization of the exciton CM wave 
function in k|| space, as the QD size increases [Eq. (|^ ]. 
This effect is consistent with the strong increase in QD 
exciton lifetime as compared with QW excitons. The 
competition between the cavity and localization effects is 
shown in the figure by the dip in the calculated emission 
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FIG. 3: (color online) Emission rates into a collection angle 
13.2° vs. exciton energy. Full (red) circles are the measured 
values given in counts per second (right axis). Solid (blue) line 
shows the calculation where detuning is obtained by varying 
the QD radius (given in the top axis). Dashed (green) line 
shows the calculation with QD size fixed to _Rd = 200A and 
the detuning obtained by varying the QD material bandgap. 
A scaling factor of 3000 was taken between the left and right 
axes (see text). 



rate (solid line) for i?D ~ 300 A (negative detuning). For 
yet larger QDs the localization effect takes over and the 
emission rate rapidly increases (the figure includes QDs 
with i?D < 6OOA) . Such large QDs are less likely to form 
in our case and this increase is therefore not observed. 

Figure^lshows the dependence of the QD spontaneous 
emission on the excitation power. At low excitation the 
PL is dominated by a sharp emission line arising from 
the recombination of the Is exciton. As the excitation 
power is increased, the dot populates with more carriers 
and recombinations of multiexcitons from higher collec- 
tive states appear in the spectrumSLSS. For still higher 
excitation powers the emission is broadened and is ex- 
tended beyond the cavity mode linewidth (see the reflec- 
tivity curve in Fig. ^ . This is qualitatively confirmed by 
calculated emission rates that are superimposed over the 
high-excitation-power PL curve (topmost PL spectrum 
in Fig. 01). As before, the detuning of higher dot states 
is represented by either varying the dot size (solid blue 
line) or its material bandgap (dashed green line). 

In Figure |21 we plot the experimental and calculated 
exciton lifetimes as functions of the emission spectral en- 
ergy. Each point in the figure represents a measurement 
from a specific QD (open circles). The sample temper- 
ature during the measurements was 15K. The emission 
lifetime was estimated from a single exponential fit to 
the initial part of the measured PL decay curves (not 
shown). In general, the decay curves were quite sensi- 
tive to the excitation intensity and they were not single- 
exponentials. For the measurements we used excitation 
intensities which exactly saturate the PL emission. We 
note that there is a considerable scatter in the measured 
data. We believe that this scatter is mainly due to non- 
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most 7SE ^ 1.6^eV ^ 7c thus we are always in the weak 
coupHng regime where our approximate solution is valid. 
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Energy (eV) 

FIG. 4: (color online) Measured reflectivity (upper line) and 
integrated PL intensity for various excitation powers. Calcu- 
lated emission rates with detuning by varying the QD radius 
(solid blue line) and its material bandgap (dashed green line) 
are superimposed on the measured PL with the highest exci- 
tation power. 



radiative processes, which may depend on the particular 
environment and charge state of a given QD. Also, as 
mentioned above, the measured decay times were excita- 
tion intensity dependent. Our particular choice of exci- 
tation intensity, for which the exciton PL saturates, may 
be somewhat arbitrary. The saturation excitation may 
vary from one QD to the other in a way that has little to 
do with the microcavity. 

For the model in which the detuning is represented by 
varying the QD material bandgap (dashed line), above 
resonance there are always cavity modes in the 2D con- 
tinuum that can couple to the exciton due to the photon 
in-plane dispersion; thus the coupling remains effective 
even for a large positive detuning, resulting in a short re- 
combination time. The emission from negatively detuned 
QDs is largely inhibited consistent with the experiment 
where emission from these QDs is not observed. This 
inhibition depends strongly on the cavity Q-factor and 
would be much less pronounced for weaker cavities. The 
calculations in which detunings are represented by vary- 
ing the dot size (solid line) show a different behavior. 
When the QD radius is reduced (energies above reso- 
nance in Fig. O the lifetime increases due to the spread 
of the exciton CM wavefunction in k\\ space, as explained 
above. The very different lifetime behavior of the two 
detuning mechanisms may aid in identifying the various 
QDs. The results of both figures 13 and [3 seem to indicate 
that the majority of QDs have similar composition and 
their variations are mainly in size. 

We note that the increase in radiative lifetime with re- 
duction of the QD size is suppressed in the limit of strong 
confinement. This is because the oscillator strength (/||) 
increases due to the localization of the relative motion 
wavefunction. This effect is superseded by the break- 
down of the selection rule caused by the spread of the 
CM k-space wavefunction outside gx- Also note that the 
calculated linewidths due to spontaneous emission are at 
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FIG. 5: Exciton radiative lifetime vs. its energy: open (red) 
circles are the measured values. Solid (blue) line shows the 
calculation where detuning is obtained with variable QD ra- 
dius (given in the top axis). Dashed (green) line shows the 
calculation with QD size fixed to Ru = 200 A and its bandgap 
varied. 

Next, angle resolved spectroscopy was performed on a 
single QD using high power excitation (8/iiy). In these 
conditions the QD is populated with many carriers re- 
sulting in many available recombination channels. The 
PL spectrum therefore reflects the density of electromag- 
netic modes in the cavity. Fig. |^ shows a contour plot 
giving the measured PL intensity of a highly excited sin- 
gle QD, as registered by the CCD array camera at the 
exit of the monochromator. The horizontal axis gives the 
spectral dispersion while the vertical axis gives the angu- 
lar distribution. Since the kz part of the emission is fixed 
by the resonance condition to the cavity mode, this angle 
corresponds to the in-plane wave vector of the emission. 
From the geometry and lenses used for the imaging we es- 
timate that each row in the CCD camera is equivalent to 
0.38mrad. The figure clearly demonstrates the in-plane 
dispersion of the cavity mode. The peaks in the emission 
intensity shown in the figure are most likely related to re- 
combinations from various carrier configurations, e.g. the 
neutral and charged excitons, biexcitons and higher order 
multiexcitons ^i*^, of the single QD. These multiexciton 
lines appear due to the high excitation power used in this 
experiment. Since our model accounts only for the neu- 
tral exciton line, we cannot compare our experimental 
data directly with the model calculations. Nevertheless, 
by taking the derivative of Eq. (|26|l with respect to the 
k|| space area (using q — q^sinO), we can achieve a sim- 
ilar in-plane dispersion of the cavity mode, as shown in 
Fig. [71 The models of QDs in which the size and the 
material bandgap were varied to represent the detuning 
show similar dispersion behavior. The difference between 
the two cases is mainly that in the case of fixed QD size 
the emission extends over a larger energy range (Fig.Et*), 
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which can be explained using the preceding arguments. 
We attribute the difference between the symmetric shape 
of the calculated angle resolved spectrum and the asym- 
metric shape of the measured one to the finite aperture 
of the confocal setup. The data was taken from a sin- 
gle QD that was positively detuned with respect to the 
MC resonance. Due to the relatively large spatial distri- 
bution of the emission from the QD (see section lIVBp , 
we could not collect the light from both negative and 
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FIG. 6: (color online) Contour plot of the PL intensity from 
a highly excited single resonant QD in a planar microcavity 
as a function of photon wavelength and emission angle. Each 
contour represents 10% of the maximal intensity. 
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FIG. 7: Contour plot of the calculated emission rate distribu- 
tions as a function of photon wavelength and emission angle 
(k||). Each contour represents 10% of the maximal rate. De- 
tuning is obtained by (a) varying the QD radius, and (b) 
varying the QD material bandgap and fixing the dot radius 
to _Rd = 200l. 



B. Spatial distribution of the Emission 

We now turn to examine the in-plane spatial spread 
of a single QD emission. Fig. shows a spatially in- 
tegrated emission spectrum from a particular QD under 
moderate excitation power. Few discrete spectral lines 
are observed. Each spectral line corresponds to a par- 
ticular emission line arising from the recombination of a 
ground state e-h pair from different many carriers' col- 
lective statesi^ For our purpose it suffices to note that 
each line is spectrally detuned differently from the cavity 
mode. In order to estimate the spectral detuning of each 
line the reflectivity spectrum from this sample is overlaid 
on the emission spectrum in Fig.|SJi. The least and most 
detuned lines are marked in the figure by A and by B 
respectively. Wavelength selective spatial images of the 
least (A) and most (B) detuned spectral lines are shown 
in Fig. and Fig. IS]:, respectively. The intensity distri- 
butions of the electromagnetic fields associated with the 
emission lines are essentially symmetric, though some ob- 
scuration, which can be partially attributed to mechan- 
ical drift, is present. The emission patterns contain a 
central strong spot surrounded by concentric rings with 
decreasing intensities. The intensities along the diago- 



nals, which are marked by dashed (red) and solid (blue) 
lines in Figs. IHb and |Ht, respectively, are displayed as 
functions of the in-plane distance from the image center 
in Fig.ili. 

Figs. EJi and |^ show the calculated images for the 
detunings of line A and line B respectively. Since both 
spectral lines originate from the same QD, we used differ- 
ent bandgaps to obtain the detuning. The images were 
generated by calculating the in-plane Fourier transform 
of the integrand in Eq. [23 The rates were calibrated so 
that their spatially integrated emission would match the 
calculated values given in Fig. |31 Fig. ^ shows the in- 
tensities along the diagonals of the calculated images as 
functions of the in-plane distance from the image center. 
In both measured and calculated images it is evident that 
when the QD is close to resonance with the cavity mode, 
its emission exhibits an in-plane spatial spread that is 
much larger than the QD size. The closer the emission 
energy to resonance is the larger is its spatial extent, and 
vise versa, the larger the detuning is, the more localized 
the emission pattern becomes. Calculated images, where 
detuning was modeled by varying the dot size, yielded 
similar behavior of the emission pattern, implying that 
the PL spatial spread is mainly a cavity effect depending 
only on detuning. In general, the calculated images agree 
with the experimental observations, both in their spatial 
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FIG. 8: (color online) (a) Spatially integrated emission spec- 
trum measured from a single QD under moderate excitation 
power. The dashed (green) line shows the normalized reflec- 
tion of the sample; Wavelength selective spatial image mea- 
sured for (b) spectral line A; and (c) spectral line B; (d) 
Emission intensities along the diagonals which are marked by 
dashed (red) and solid (blue) lines as functions of the distance 
from the image center. 
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FIG. 9; (color online) Calculated spatial spread of the emis- 
sion from a single QD detuned from resonance by (a) S — 
3meV; (b) 5 — 7.4meV, corresponding to the measured spec- 
tral lines marked A and B in Fig. |H| respectively. Detuning 
is modeled by varying the QD material bandgap while the 
QD radius is fixed to Ru = 200 A); (c) Calculated emission 
intensities as functions of the distance from the image center. 
(Constant background of 400 counts was added to the cal- 
culated emission to facilitate comparison with the measured 
data). 



extent and in their dependence on the detuning, as can 
be concluded by comparing Fig. |H1 with Fig.|Sl 

A discernible feature in Figs. (SJi andl^t is the depen- 
dence of the spacing between the rings in the PL images 
on the detuning. The rings become more pronounced 
and denser for increasing positive detuning, where larger 
k\\ are needed for the emission. This gives rise to oscilla- 
tions that are superimposed on the central emission line 
due to contributions arising from larger fc|| modes. The 
dependencies of the calculated locations of the 2nd and 
3rd peaks on detuning are shown in Figure lTUk . and their 
intensities normalized to the central peak are shown in 
Fig.llOb. confirming these observations. Identical results 
are obtained for both detuning models, suggesting again 
that this emission pattern is a cavity effect, depending 
only on detuning. Together with the central peak spread, 
the locations of the emission rings may aid in evaluating 
the QD's detuning. We note, though, that as Fig. IlUb 
suggests, for spectral lines that are close to resonance 
with the cavity mode, the relative intensity of the rings 
emission to the central spot decreases rapidly, making 
them difficult to measure. 
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FIG. 10: (color online) (a) Calculated positions of second 
and third rings in the emission distribution; (b) Calculated 
emission intensity of second and third rings normalized to 
central emission peak. 
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V. SUMMARY 

In this work we have studied the interaction between 

the electromagnetic modes in a planar cavity and cx- 
citons in single QDs embedded at the antinode of the 
cavity. 

We find that the spontaneous emission outside the mi- 
crocavity is enhanced considerably for dots in resonance 
with the cavity mode. The emission intensity rapidly 
decreases with detuning for lines that are positively de- 
tuned from the cavity mode. Emission from negatively 
detuned QDs is largely inhibited, and we were unable to 
observe any spectral emission below the cavity mode en- 
ergy. Lifetime measurements have shown a decrease by a 
factor of roughly two for the lifetime of resonant QDs as 
compared with QDs that are 20 meV detuned from the 
cavity line. In-plane dispersion of cavity photons were di- 
rectly observed by strong excitation of the QD, thereby 
transforming it into a broadband light source. Finally, 
we obtained near field images of the spatial emission dis- 
tributions from single QDs, showing a large spread over 
several microns in the case of resonant QDs. These im- 
ages have also shown emission rings with separation and 



intensity that depend on the detuning of the QD from 
the cavity mode. 

All of these features were qualitatively accounted for 
by calculations where detuning from the cavity mode was 
modeled by either varying the QD size or its bandgap. 
For most of the experimental data we have found bet- 
ter agreement with the first mechanism, suggesting that 
most of the QDs in strain driven self assembled sam- 
ples vary in size rather than in composition. The range 
of 90A < Rd < 300A that accounts for the observed 
phenomena is consistent with typical self assembled In- 
GaAs/GaAs quantum dots. 
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Appendix 



Here we provide some details of the calculation of the exciton's relative motion energy. 



Er = {XeXh(t>\Hr\(t>XhXe 



= 27rAC' 



^^^2a-lg-2(r/r,)« + ^2^2 



L Vh 



47r(A/;A4eA4J' / drdqrMqr)e-'^'/^^'' Y.Ul) + (Y^^ j'^^^ 
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fi = 4.k1 + q^ + q^ cos(fcjL) , gi = 2qki sm{kiL). 
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